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Abstract: We solve the loop equations of the hermitian 2-matrix model to all 
orders in the topological 1/N 2 expansion, i.e. we obtain all non-mixed correlation 
functions, in terms of residues on an algebraic curve. We give two representations of 
those residues as Feynman-like graphs, one of them involving only cubic vertices. 
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1. Introduction 

The purpose of this article is to generalize the method invented in [13], for the 
2-matrix model. The method of [13] is a diagrammatic technique for computing 
correlation functions of the 1-matrix model in terms of residues on some algebraic 
curve. 

Random matrix models play an important role in physics and mathematics [32], 
and have a wealth of applications which are too long to list here. In this article, 
we consider "formal" random matrix integrals, which are known to be generating 
functions for counting some classes of discrete surfaces [7,10,21,22,34]. 

The partition function, free energy and correlation functions are all generating 
functions enumerating some kinds of graphs (respectively closed graphs, connected 
closed graphs, open graphs), which graphs can be seen as discrete surfaces. 

In the formal model, the size N of matrices, is just a complex parameter, it needs 
not be an integer, and all observables (free energy, correlation functions) always have 
a 1/N expansion, because for each power of the expansion parameters, there is only 
a finite number of graphs with a given power of N. The power of N in a graph is 
its Euler characteristic, and thus the 1/N expansion is known as the "topological 
expansion" discovered by 't Hooft [34]. In the formal model, N is thus an expansion 
parameter, and working order by order in N enumerates only discrete surfaces of 
a given topology [7]. An efficient method for dealing with this formal model is to 
consider the Schwinger-Dyson equations, called loop equations in this context [10,33]. 

To large N limit (i.e. planar topologies), the solution of loop equations is known 
to be related to Toda hierarchy [9,28,35,36]. For this reason, the large N expansion 
of matrix models plays an important role in integrable systems, and in many areas of 
physics [29]. It was understood by [11] that the low energy effective action of some 
string theory models is also described by matrix models. 
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In the beginning, formal matrix models were considered omlyin their 1-cut phase, 
because a potential which is a small deformation of a quadratic one, must have only 
one well, i.e. the variables perturbatively explore only one well. However, a N x N 
matrix has N eigenvalues, and even though each of them can explore perturbatively 
only one well, they do not need explore all the same well. That gives "multicut" 
solutions of matrix models, where the number of eigenvalues near each extremum of 
the potential is fixed (fixed filling fractions). Multicut solutions play an important 
role in string theory, as they describe multi-particle states [11,12]. Multicut solutions 
correspond to enumerating surfaces with contact terms, which can be called "foam 
of surfaces" as described in [6,22]. 

The link between formal matrix models (which always have a 1/JV expansion) 
and convergent matrix integrals (which have a 1/N expansion only in the 1-cut case 
under certain assumptions), has been better understood after the work of [6]. We 
emphasize again, that the results developed in this article concern the formal matrix 
model with fixed filling fractions, and should not be applied to convergent matrix 
model directly. 

Recently, it has progressively become clear that large N expansion of random 
matrix models has a strong link with algebraic geometry [27]. The free energy and 
correlation functions have been computed in terms of properties of an algebraic 
curve. The large N limit of the 1-point correlation function (called the resolvent) is 
solution of an algebraic equation, which thus defines an algebraic curve. There have 
been many works which computed free energy and correlation functions in terms of 
that algebraic curve. The leading order resolvent and free energy were computed in 
the 1-cut case (algebraic curve of genus zero) in the pioneering work of [7], then some 
recursive method for computing correlation functions and free energy to all orders in 
1/N were invented by [3,4]. Those methods were first limited to 1-matrix case and 
1-cut. 

Then for 1-matrix several works have dealt with multicut: Akeman and Ambj0rn 
found the first subleading term for the multicut resolvent and the 2-cut free energy 
[1,2], Chekhov [8] and one of the authors together with Kokotov and Korotkin [17] 
found simultaneously the first subleading term for the multi-cut free energy . Then 
a (non-recursive) diagrammatic method was invented in [13] to find all correlation 
functions to all orders, in the multicut case. 

The 1-matrix model, corresponds to hyper elliptical curves only. In order to have 
more general algebraic curves, one needs at least a 2-matrix model. For the 2-matrix 
models, the loop equations have been known since [33], and have been written in 
a concise form in [18-20]. They have been used to find the subleading term of the 
free energy, first in the genus zero case in [15], then in the genus 1 case in [16], 
and with arbitrary genus in [17]. The purpose of this article is to generalize the 
diagrammatic method of [13] for the computation of non-mixed correlation functions 
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in the 2-matrix case. We solve the loop equations and present their solutions (the 
non-mixed correlation function's expansion) under two different diagrammatic forms. 
We first build a cubic diagrammatic representation before presenting an effective non 
cubic theory. 

Outline of the article: 

• In section 2, we introduce the model and our notations. 

• Section 3 is dedicated to the derivation of loop equations. We derive the funda- 
mental "master loop equation" before deriving loop equations whose solutions 
are non-mixed correlation functions 

• In section 4, we show how a compact Riemann surface arises from the leading 
order of the master loop equation and present notations and tools of algebraic 
geometry needed for the computation of correlation functions. 

• In section 5, we present a diagrammatic solution of the loop equations as cubic 
Feynman-like graphs. 

• Section 6 is dedicated to the presentation of another representation of the non- 
mixed correlation functions as graphs of a non cubic effective theory. 

• In section 7, we study the example of the gaussian case corresponding to the 
1-matrix model limit. 

2. Definitions and notations 

2.1 Definition of the formal 2-matrix model with fixed filling fractions 

In this article, we are interested in the study of the formal-two-matrix-model and 
the computation of a whole family of observables. The partition function Z is the 
formal matrix integral: 

Z:=f _ dMldM2e -NTr(y,{M l)+ V 2 {M 2 )-M,M 2 ) ^ 

where M\ and M 2 are two N x N hermitian matrices, dM\ and dM 2 the products of 
Lebesgue measures of the real components of M 1 and M 2 respectively, and V\ and 
V 2 two polynomial potentials of degree di + 1 and d 2 + 1 respectively : 

di+l d 2 +l ~ 

k=i K k=i K 

Formal integral means it is computed as the formal power series expansion order 
by order in the g^s (see [7,10,34]) of a matrix integral, where the non-quadratic terms 
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in the potentials V\ and V 2 are treated as perturbations near quadratic potentials. 
Such a perturbative expansion can be performed only near local extrema of V\(x) + 
Viiy) — %U, i-e. near points such that: 

Vm=m , Vi(rH)=ti (2.3) 

which has did 2 solutions. Therefore, if M 1 and M 2 are diagonal matrices, whose di- 
agonal entries are some £j's (resp. 77* 's), (Mi, M 2 ) is a local extremum of tr(Vi(Mi) + 
V2(M 2 ) — M1M2) around which we can perform a perturbative expansion. 

The choice of such an extremum, around which the perturbative series is com- 
puted, is equivalent to the choice of the number of eigenvalues near each pair (£3, rji), 
i — 1, . . . , did 2 , i.e. the data of d\d 2 integers rii such that: 

E n t = N (2.4) 

i=i 

This means, that we can choose some contours C iy % = 1, . . . , did 2 , such that the 
following equality holds order by order in the perturbative expansion: 

The numbers ^ are called filling fractions. Thus, in the formal model, filling 
fractions are fixed parameters. 

Fat graphs and discrete random surfaces 

Once filling fractions are chosen, we perform the perturbative expansion. Each 
term of that formal expansion is an expectation value of a gaussian integral, and 
using Wick's theorem, each term can be represented by a Feynman graph. Because 
the integration variables are matrices, the graphs are "fat graphs", which have a 
2-dimensional structure. The Hermitean matrix models thus enumerate oriented 
surfaces (other matrix ensembles can enumerate non-oriented surfaces). This Formal 
expansion equivalent to an enumerating function of Feynman graphs is a standard 
tool in physics [10,34]. Random matrices have thus played a role in all theories 
where one needs to sum over surfaces, i.e. string theory and quantum gravity (i.e. 
statistical physics on a random lattice). 

Following this interpretation, the loop equations [33] can be understood as rela- 
tionships linking surfaces of different genus and different number of boundaries. 

2.2 Notations 

2.2.1 Notation for sets of variables 

We will consider functions of many variables x±,x 2 , x^, . . . , Xk, or of a subset of those 
variables. In that purpose we introduce the following notations: 
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Let K be a k— upple of integers: 

K = (i 1 ,i 2 ,---,ik) (2.6) 

We denote fc = the length (or cardinal) of K. For any j < \K\, we denote Kj 
the set of all j— upples (i.e. subsets of length j) contained in K: 

K r .= {JcK , \J\=j} (2.7) 

We define the following A;— upple of complex numbers: 

Xjf . i j • • • j ) (2-8) 

2.2.2 Correlation functions 

For a given fc, we define the correlation function: 

w k { Xl , ...,x k ):= N k ~ 2 (il tr 1 \ , (2.9) 

\i=l i/c 

i.e., with the previous notations: 



^,(x^):=iV^- 2 (ntr-4^) , (2-10) 

\r=l x *r JW 1/^ 



where the formal average (.) is computed with the measure in eq. fl2.1|) , and the 
subscript c means connected part (cumulant). 

Those correlation functions can be expanded as formal series in in the large 
N limit: 

«J fc (x A -) = E^4 ft) (xi.) (2.11) 

The purpose of this article is to compute wf'(x;c) as residues on an algebraic 
curve and represent it with Feynman-like graphs of a cubic field theory on the curve. 

We also define the following auxiliary functions: 

Mx,V,*k) == Wtr \, Vi(v) ~ f M2) S tr — L_\ (2.12) 

\ x-Mi y-M 2 jL\ x lr -hhj x ' 



x-Mx y-M 2 ^ x ir -M 1 / c 

(2-13) 



5 



Notice that u k ,( x i V't is a polynomial in y of degree d 2 — 1, and p k (x, y\ x#) is a 
polynomial in x of degree d\ — 1 and in y of degree d 2 — 1- 
It is convenient to renormalize those functions, and define: 



u k (x, y; x. K ) := u k (x, y; x K ) - 5 kfi (V 2 (y) - x) (2.15) 

and 

w k (yi K ) := W k (yL K ) + 7 ' — (2.16) 

(xi - x 2 y 

Let us remark that all those functions have the same kind of topological ex- 
pansion as Wfc(x^) and one defines p k h \x,y;xx) and u k h \x,y;^K) as well like in 
eq. (|2XID. 



We define the function: 

Y(x) := V[(x) - Wl (x) (2.17) 

which we see below, describes the algebraic curve. 

The -^2 expansion of such correlation functions is known to enumerate discrete 
surfaces of a given topology, whose polygons carry a spin + or - (Ising model on a 
random surface [26,29]), see [22] for the multicut case i.e. foam of Ising surfaces. 

The w k are generating functions enumerating genus h discrete surfaces with k 
boundaries of spin +. 

As an example, w 2 enumerates surfaces of genus 3 with 2 boundaries: 



_(3) 




(2.18) 



Notice that the question of boundaries with non uniform spin, i.e. with changes 
of boundary conditions has been solved to leading order only in [14]. 

3. Loop equations 

There exist several methods for computing the free energy and correlation func- 
tions, the one we consider here is the "loop equation" method, which is nothing 
but Schwinger-Dyson, or Ward identities [10,33]. They implement the Virasoro or 
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W-algebra constraints on the partition function [27,31], i.e. the fact that the matrix 
integral is left unchanged under a change of variable. The loop equations are valid 
in the formal model, order by order in the expansion parameters. 

For the 2- matrix model, loop equations have been known since [33] , and written 
in a more systematic way in [18-20,27]. 

3.1 The master loop equation 

It is well known that in the large N limit, loop equations imply an algebraic equation 
for the functions w±, i.e. for the function Y(x), called the master loop equation. Let 
us briefly recall how to derive it (see [20]): 

• the change of variables M 2 — > M 2 + £ X 1 M implies: 

= a (x) — xw\{x) + 1 (3.1) 

• the change of variables Mi — > Mi + e x 2 Ml ^^zi|~^ implies: 

w 1 (x)u Q (x, y) + V, x) = V((x)u (x, y) - p (x, y) - yu (x, y) 

+V^{y)w 1 {x) - a (x) 

(3.2) 

i.e., putting everything together: 
(y - Y(x))u (x,y) + ^Mi(x,y;x) = (V£(y) - x)(V{(x) - y) -p (x,y) + l (3.3) 
We define: 

E(x,y) = (Vl(y) - x)(V((x) - y)- Po (x,y) + 1 (3.4) 
The master loop equation is thus: 



(y - Y(x))u (x, y) + -^u^x, y- x) = E(x, y) 



(3.5) 



where E(x, y) is a polynomial of degree di + 1 in x and d 2 + 1 in y. 
3.2 Loop equations for correlation functions 

We now derive the loop equations which allow to compute recursively the k-point 
non-mixed correlation functions. 

• The change of variables 5M 2 = x } Ml Vu=\ tr x .l Ml implies (see [20]): 

a k (x\ x x ) = x w k+1 (x, xjf) - N 2 w k (y. K ) (3.6) 
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The change of variables 8M\ 



[20]): 



XlU tr ^ implies (see 



fc-1 

u7i(x)M fc (a;,?/;x K ) + ^ Uj(x, y- Xj) HJ fc _ i+1 (x, x#_j) 
j=o Jei^j 

1 - / 

+ \p gfc-i(x, xa-,^-}) - Hk^jxj, y; * K -{j}) 

j=i x x j 

VjC^Wfe.oC^, 2/; xj?) - p k (x, y; x#) 
-yu k (x, y\ x K ) + V% (y)w k+ i(x, x K ) - afc(s; xjc) 



(3.7) 
i.e. for k > 1: 



(y - ^(^))«fc(^,2/;xi<r) 



- Ej=o Ejex, ujfo y; xj) w fc _ i+ i(x, x^_j) 
-^M fc+ i(x,y;x,x x ) 



(3.8) 

The purpose of this article is to solve eq. (|3.8|) and compute w k for all k and h. 

4. Leading order and algebraic geometry 

4.1 Leading order of the master loop equation 

To large N leading order, the master loop equation eq. ( |3.5| ) reads: 



(y - Y(x))u (x, y) = E(x, y) 

(4.1) 

Since uq(x, y) is a polynomial in y, it has no singularity for y finite and the LHS 
vanishes for y = Y(x), i.e.: 



E(x,Y(x)) = 

This defines an algebraic curve E(x, y) = 0. 
Notice that to leading order we have: 

E(x,y) 



(4.2) 



u {x,y) 



and 



y ~ Y(x) 
u (x,Y(x)) = E y (x,Y(x)) 



(4.3) 
(4.4) 
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4.2 Introduction to some algebraic geometry 

We use notations similar to [25] or [24]. Some useful hints for understanding this 
section can be found in Appendix A. 

Let us parameterize the curve E(x,y) = with a running point p of a compact 
Riemann surface £. It means that we define two meromorphic functions x(p) and 
y(j)) on £ such that: 

E(x, y) = 3p e £ x = x{p) , y = y{p) (4.5) 

The functions x and y are not bijective. Indeed, since E(x,y) is a polynomial of 
degree d 2 + 1 in y, it has d 2 + 1 solutions, i.e. for a given x, there exist rf 2 + 1 points 
p on £ such that x{p) = x. Thus, the Riemann surface is made of d 2 + 1 ^-sheets, 
respectively d\ + 1 y-sheets. Hence, from now on, we use these notations: 

x(p) = x p = p 3 (x) for j = 0,...,d 2 (4.6) 

y(p) =y <&p=p>(x) for j = 0, . . . ,di (4.7) 

We will most often omit the exponent corresponding to the physical sheet: 
p = p°. 

For instance, one can write E(x,y) as: 

di 

E(x(p),y(q)) = -g dl+ i x JJ(a;(p) -x(q l (y))) 

j=0 

-faixn^-^W)) (4-8) 

i=0 



Considering that the tuj^'s, uj^'s and p^'s are multivalued functions in their 
arguments x, we now work with differentials monovalued on the Riemann surface. 
Let us write the differentials: 



W k+1 {p, p K ) := w k+1 (x(p),x(p K ))dx(p) J] dx (Pi) ( 4 - 9 ) 

i=i 

k 

Uk(p, V, Pk) ■= u k (x(p),y; x(p K ))dx(p) dxfa) (4.10) 

i=i 

fe 

P fe (x,y;p K ) :=p fe (x,y;x(p K )) JJ tte(pi) (4.11) 

i=i 

Note: In the following, the arguments of a function will be called x(p) or y(r) 
if the function is defined on the basis, and p or r if the function is defined on the 
Riemann surface - and so multivalued on the basis-. 
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Let us now review the notations we use in this article to denote some basic 
objects. For definitions and details, we refer the reader to Appendix A and [25] 
or [24]. 

•Canonical cycles: Ai, B{ for i = 1, . . . , g where g is the genus of the compact 
Riemann surface S (0 < g < d\d 2 — 1), such that: 

AiHBi = 6i d (4.12) 

•Branch points in x: They are the zeroes of dx on the surface. We denote 
them by ai, i = 1, . . . ,d 2 + 1 + 2g. 

•Bergmann kernel: It is the unique bilinear differential with only one double 
pole at p = q satisfying: 

B(p,q) ~ d <^ dx ^ + finite and y t I B(p,q) = (4.13) 
p—9 (x(p) - x{q)Y JptAi 

•Abelian differential of third kind: It is the differential defined by dS q ^ r (p) = 
J ql=r B(p, q'). Notice that it has the following properties: 

Res dS q>r (p) = 1 = - Res dS q , r {p) and Vz jf dS q>r (p) = (4.14) 

4.3 Fixed filling fractions 

To large N leading order, the loop equation eq. Q4.2| ) is an algebraic equation: 

E(x,Y(x)) = (4.15) 

The coefficients of E are determined using filling fractions. Since Wi(x) = V{(x) — 
Y(x), eq. (|2.5|) gives (up to a redefinition of Ai): 

/ ydx = / xdy = e« (4.16) 

2m JAi 2?7r J Ai 



Let us recall that (see section [2.1|) the e^'s are called filling fractions, and they 
are given parameters (moduli) of the model. They don't depend on the potential or 
on any other parameter. 

In particular, since all correlation functions w^{x\, . . . , xj.) are obtained by deriva- 
tion of W\ with respect to the potential V\ ( [4]), we have for k > 2: 

Wk(xx, ...,x k )dxi = (4.17) 



2m JAi 

Equation eq. (|4.16|) together with the large x and y behaviors eq. (|A.2p and 
eq. ( |A.1|) , are sufficient to determine completely all the coefficients of the polynomial 
E(x,y), and thus the leading large N resolvent Wi(x). 

In what follows, we assume that the leading resolvent, i.e. the function Y(x) is 
known, and we refer the reader to the existing literature on that topic, for instance 
[5,20,27,30]. 
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5. Diagrammatic solution as cubic graphs 

In this section we present a first way of describing the solution of the loop equation 
eq. ( |3.8| ) by trivalent diagrams whose h loop level corresponds to the h-th term 
of the topological expansion. 

5.1 Solution in the planar limit 

Before considering the full expansion, let us focus on the structure of the leading 
terms corresponding to planar fat graphs. Thus the 1/N 2 terms in the loop equations 
are omited. 

From now on and particularly in this paragraph, we drop the genus zero exponent 
(0) when it is clear that we deal with the planar limit, i.e. w^\x.k) — > u>fc(x#). 

Up to now, the loop equations were written in terms of multivalued functions. 
It is more appropriate to write them in terms of meromorphic differentials on the 
Riemann surface. Thus, one writes eq. ( |3.8| ) in the planar limit as follows: 



(y(r)-y(p))U k (p,y(r);p K ) = -EtSEje*. 



Hp Uj(p,y(r);pj)W k - ]+1 (p,p K 



3=0 ^JeKj dx ( p ) 

, V fc j ( Uk-i{Pj,y{r);v K -{ 3 }) dx( P ) \ /c-i\ 
^2^ j=1 a Pj y x( P )-x{ P] ) dx( Pj )J ^- ) 

-P k (x(p),y(r);p K )dx(p) 

Starting from eq. ( |5.1|) , we determine W k and U k for any k by recursion on k. 

Let us assume that one knows Wj(pj) for j < k and Uj(p, pj) for j < k — 1. The 
first step consists in the determination of Wk+i(p, Pk) as a function of the lower order 
correlation functions. The second step leads to the computation of U k (p, Pk)- Once 
this is done, one knows the correlation functions one order upper. The initial terms 
W2 and U\ can be found in the literature [5,20,27] and are rederived in Appendix B. 

5.1.1 Determination of W k +i for k > 2 

If one chooses r = p in eq. (|5.1|) , one gets (using eq. ( |4.3|) and eq. 



Ey{x(p),y{p))W k +i(p,PK) = -P k (x(p),y(p);p K )dx(p) 

Uj(p,y(p);pj) ^k- 3 +i{p,PK-j) /k q\ 

,spk j f Uk-i(pj,y(p);PK-{j}) dx( P ) 
-r 2^j=i u P j y x {p)-x{ Pj ) dx( Pj ) 

Notice that the two equations eq. (|5.1| ) and eq. (|5.2|) imply by recursion, that 
W k and Uk are indeed meromorphic differentials on the curve, in all their variables. 



We define: 



fc E y {x{p>),y{j?))dx{p>) v ; 
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Note that we have already obtained (see eq. ( |4.3| )) that: 

Ri(p l ) = 



Using eq. (|4.14j ), the Cauchy formula gives: 

W k+1 (p,p K ) = - Res Wk+i(p' ,-pK)dS p r 7 o(p) 



(5.4) 



(5.5) 



where o G £ is an arbitrary point on the Riemann surface. 

The integrand has poles in p' only at p' = p and the branch points p' = a s (this 
can be proven recursively by differentiating wrt the potential ^fr)- Using Riemann 
bilinear identity eq. ( |A.10| ), we can then move the integration contour and get: 

W k+1 (p,p K ) = V Res W k+1 (p',p K )dS p > >0 (p) (5.6) 

We now introduce the loop equation eq. ( |5.2| ) inside this expression and remark 
that only one term has poles when p' — > a s . Thus W k +i(p, Pk) can be written: 



W k+1 {p, p K ) 



£ Res £ E U /^}^ 



k-l 



E y (x(p'),y(p')) dx(p > ) 



E Res E E R°(p',Pj)w k - j+1 (p', PK „j)ds pl M 



p'^a 



(5.7) 



3=1 J & K 3 



Notice that U k (p,y;pK) is a polynomial in y whose degree is equal to di — 1. 
Considering its di values for y = y{p l ) with i G [l,^], the interpolation formula 
reads: 



Wy 



{y -y(p))u k {p,y,p K ) 



E(x(p),y) 
for y = y(p), this gives: 



Y u k(Piy(p%PK)(y(j>) -yjf)) 
h {y-y{p i ))E y {x{p),y{p i )) 



<h 



R° k (p,p K ) = -E^fop*) 



(5.9) 



i=l 



So, in eq. ( |5.7|) , one obtains the recursive formula for W k (px)'- 



d 2 k-l 



w k+1 (p,p K ) = EE E E Res RW-,vjW k - j+1 (p\p K -j)ds p ,M 



i=i 3=i jeKj 



p'-^a. 



(5.10) 



The sum over j represents the summation over all partitions of K into two subsets 
J and K — J. 
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5.1.2 Determination of R l k 

In this section, we find a recursion formula for R\. 

For this purpose, one needs to know an intermediate expression defining the 
different U^s as well as a relation linking the value of 

k— 1 

E Uj(p\ y( P ); Vj)W k - j+l { P \ pa-,/) (5.11) 

3=0 

for different i's. 

Let us rewrite here eq. (|5.1j) : 



fc-i j 

(j/(r) -y(q))U k (q,y(r);p K ) = ~E E -fT\ U ^ q ' Pj "> Wk ~i+^ Pk-j) 

3=0 J&K 3 ax \Q) 

+ Y^ d f U k-i(Pj^( r )^PK-{j}) dx(g) \ 



3=1 



x(q) — x(pj) dx(pj) J 



-P k (x(q),y(r);p K )dx(q) (5.12) 

In what follows, we use the properties of rational functions defined on the basis 
and not on the Riemann surface (for some more details, see the case k — 1 in 
Appendix B). 

For r = q = p\ eq. ( |5.12j ) reads: 



k—l i 

o = -E E T^ u 3(p l MpJ,pj)w k ^ +l (p\p K .j) 

3=0 J£Kj ^{P J 

/ f/ fc _i(p i ,y(p l );p A '-{j}) dx(p*)\ 



3=1 



xip 1 ) — x(pj) dx(pj) ; 



-P k (x(p l ),y(p t )]p K )dx{p t 



k-l i 



i=o JeKj ax \P) 



+ E d w 

-P fc (:r(p),j/(p*);p*:)<&;(p) (5.13) 



where we have used that x{p) = x{p l ). 
Now, write eq. ( 5.12|) with r = p % and q = p: 



W) -y(p))Uh(p,y(j?);PK) 



fc-i ^ 



E E ^ife 2/(^)5 Pj)^-i+i(P>P^-j) 
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k 

+ E d Pj 

3=1 



'Uk-iiPj.vip^WK-ij}) dx(p)\ 



x(p) — x(pj) dx(pj) J 
-P k (x(p),y(p i );p K )dx(p) (5.14) 



and inserting eq. ( |5.13| ) we get: 



(vipi - y(p))u k (p,y(p i );p K ) 

k-l ^ 

= -E E -fr^UiiPivip^vj) w k „ j+1 (p,p K _j) 

j=0 JizRj ax {P) 
k—1 1 

+ E E j^Uj&ytfhpjWwtfi-pK-j) (5.15) 

j=0 JeKj ax \P) 

This formula is in principle sufficient to compute the Life's recursively, and then, 
one can compute the -R^'s. However, what we need in order to get diagrammatic 
rules, is a closed recursion relation for the R^s themselves. In order to achieve this 
aim, we show that: 

Lemma: for any k > 1, one has: 



E{x(p),y)dx(p) ^ * 

U k {p,y;PK) = n — E E E 



n 



r=l KxVi...\JK r =K hi-h^-Hr=\ 



t=i (y - y(p [jt) )) dx{p) 

(5.16) 

where the sum over K\ U . . . U K r — K is a sum over all partitions of K into r 
subsets. 

Proof: It can be proven easily by recursive action of d/dVi, as in [4], however, in 
order to have a self-contained method, we want to derive it here only from the loop 
equations eq. (|3.8|) . 

The proof works by recursion on k. It is proven in Appendix B for k — 1. Let 
us assume that, it holds for any I < k — 1. 

Notice, that since both sides of eq. (|5.16|) are polynomials of y, of degree d 2 — 1, it 



is sufficient to prove that the equality holds for d 2 values of y, namely, it is sufficient 
to prove it for y = y(p l ), i = 1, . . . , d 2 . Therefore, one has to prove that: 

Uk(p,y(f);PK) _ E y (x( P l ),y(f)) j^ y y 
W\ Kr \+i (p* , pk t ) r fr W lKtl+1 {pi\p Kt ) 



dx (p) fJi (y-y (p* ))dx{p) 
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where only the sums in which one of the jts is equal to i contribute. 
The recursion hypothesis for j < k — 1, and any J G Kj gives: 



r-73 = E y {x{p),y(p))^ ^ ^ 

TT WWl(P*»P*) 

t =i (y(p*) - y(p jt )) dx(p) 

(5.18) 

In order to compute L^(p, ?/(p 4 ); pj), one has to keep only terms in the sum such 
that there exists a £ such that jt = i, i.e. 



C/ ' (P ^ iPj) - E v ( X (p'),y(p')) y, E £ 

W\J r \ +l {p\Y>J r ) g ^| Jt | + 1 (^, Pj , 



(y(p*) - <&(p) t =i (j/(p*) - yCp 7 ' 4 )) ^(p) 

(5.19) 



Insert that into eq. ( |5.15| ): 



^(p 4 ) -y(p))C/ fe (p,y(p i );P^) 

-W),^))EEi; E E 

W^-j+i(p,py-j) 7 , ' ; ^ , , . 11 



(l/(p*) - yip)) dx (p) t=i (l/W - S/(p*)) da:(p) 

+^(p i ),y(p i ))E E E E E 

3=0 JeKj r=l .hU...U.Jr = J j 1 ^j 2 ^...^jr^i 

r=i {y{p l ) - yip 3 ')) dx{p) 

(5.20) 

The difference between these two summation, keeps only j t 7^ 0,z, thus: 



U k (p,y(p % );PK) 

= E y {x{p%y{p i ))dx{p)Y^Y.J: E E 

J= o JeA' 3 r=l JiU...uJ r =J j 1 ^j 2 ^...^j r ^i,o 

(y(p*) - y(p)) t=i (vip*) - y(p jt )) dx (p) 

(5.21) 
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i.e. we have proven the lemma for k, for y = y(p l ), and since both sides are 
polynomials in y of degree g?2 — 1, the equality holds for all y. 



Theorem: For all k > 1, one has: 

Efei E.to 1 Ej e ^ W y(p); p^-i+M p*_jr) 
= EH Ej £ a- t^-(p, y(p); pj)W k - j+1 (p, p K -j) 1 ' J 

Proof of the theorem: Let us simply perform some basic rearrangements: 

d 2 fc-1 

EE E U j (p\y(p);p J )W k ^ j+1 (p i ,p K _j) 
i=i j=o jeJCj 

E E VWi(ApA-J%|+i(Ay(p);Px) 

A" 1 (jL=A'ii = 1 

fife ^2 

= E y {x{p),y{p))dx{p) E EE E E 

A-i [J L=K JX=1 r=l K 2 U...UK r+1 =L j 2 ^j 3 ^...^j r e[l,d 2 ]-{jl} 

W lKr+ll+1 {p,p Kr+1 ) " W^.I+iCp^pO 



(5.23) 



(1/(P) - Z/l^ 1 )) a= 2 W) ~ y(p 3a ))dx{p) 
d 2 d 2 

E y (x(p),y(p))dx(p)J2 E E 

r=l A-iU...UA- r+ i=A- j 1 ^j 2J t...^j r =l 

t~t ^|Jf n |+l(P ? ' n , PgJ W|jf r+ i|+l(P, PA r+1 ) 

fJi (y(p) - y(p ja ))dx(p) 

E Wj*r+ii+i(p> px r+1 )^|j|(p, y(p); pj) 
A- r+1 y j=k 



This identity simplifies eq. (|5.15|) which becomes now: 



(y^) -y(p))Rl(p,PK)dx( P ) 



fc-1 



(5.24) 



wWp\pa-) + E E E 



Uj{p l , y{p% pj)w k - j+1 {p l , pk~ 

E y (x(p),y(p i ))dx(p) 



One can now write down the final recursion formula for R\(p, Pk) in these terms: 



Rk(P, Pk) 



(y(p % )-y(p))d.x(p) 

i^k-isp sr R l j(p l ,pj)w k -3+i(p l ,PK-j) 
-r 2^j=i l^j&Kj 2^l^o,i (y(p*)-y(p))dx(p) 



(5.25) 
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The relations eq. (|5.7|) and eq. (|5.25| ) allow to compute recursively W k for any k. 
This solution can be represented by binary trees as it is presented in section ( |5.3| ). 

5.2 Solution for any genus 

In the previous paragraph, one has kept only the leading terms when performing the 
changes of variables to obtain the Schwinger-Dyson equations. Let us now write the 
corrective term for the same changes of variables so that we write a system of 
equations giving the whole 4^ expansion. One obtains the following loop equations 



(y(r) - y{p))U k {p, y{r); p K ) 

= -P k (x(p), y(r); p K )dx(p) - £j=o i^O, y(r); p,j)W k _ j+1 (p, p^-j) (5.26) 

l Uk+i(p,y(r);p,Pk) | v j ( Uk-i(pj,y(r);p K -{j}) dx( P ) \ 

N 2 dx(p) ~r u Pj \ x(p)-x(pj) dx( Pi ) ) 

For the following, one should remind the expression of the function Y(x(p)): 

Y(x) := V{(x) - w x {x) (5.27) 

Then, for h > 1: 

Y{h){x{p))= -^Wr (5,28) 

Consider now the expansion of this equation order by order. The genus h 
term (corresponding to the j^k term) gives: 

(y(r) -y(j>)M h) (p,y(ry,PK) -Y: h m= iY^(x(p))ui h - m) (p,y(r);PK) 
= -Pf> (x(p),y(r);p K )dx(p) 

£™=o E.to 1 4^j m) (p, y(r); Pj)Wtj?}(p, Pk-j) ( 5 ' 29 ) 



U ( k h +1 1) (P,y(r);p,p k ) v , / ui h \(j>j,y(r);PK-{j}) dx(p) \ 
dx(p) ' ^3 a Pj I x(p)-x( Pj ) dxipj) J 

When y(r) = y(p): 



E h m=1 Y^(x( P ))ut m \p,y(p);PK) 

P { k h \x( P ), y(p); p K )dx{p) + J2 h T=0 Efco 1 4^j m) (p, v(p); pj)w££J(p, pk- 

dx(p) Pi \ x(p)—x( P j) dx(pj) 

(5.30) 

These two equations are the generalization of eq. (5TT) and eq. ( |5.2| ) for any genus 



in the topological expansion. With all these tools, we are now able to compute all 
the terms of the expansion of non mixed traces. 

In this section, we proceed in two steps to compute the correlation function wjf 1 ^ 
for any k and any h, and represent it as a Feynman graph with h loops. The first step 
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consists in the determination of a recursive relation for Wj: , whereas the second one 

gives := J J , k , %. ^7^77^ considered the lower order terms known. 

For the following, let h and k be two given positive integers. Let us consider 
Wj™^ known for any j if m < h and any j < k if m = h. One also assume that 
^>*'( m ) j s k nown f or a ny i and any j if m < h and any j < k if m — h. Starting from 
these assumptions, one computes and R% \ what will allow to know any term 

recursively. 

5.2.1 A recursive formula for 

Let us remind eq. ( |5.30| ) in a more suitable way to emphasize that it allows us to 
compute Wfc+ifapic) with our assumption: 



W^ 1 (p,p K )U (p,y(p)) = 

-Pi h \p,y(p);p K )dx(p) 2 ( 531 ) 
- E™= u i m) (P> Pj) W k-~jTl (P. P*-j) 

-c/fc'W);^) + T., E, d P , C , '-'Sff^-°" ^i) *« 

Remark that the RHS contains only known terms except P^. h \p,y{p)',PK)- For- 
tunately, it plays no role in Cauchy formula. 

Indeed, we write the Cauchy formula, move the integration contour and vanish 
integrals around the cycles thanks to the Riemann bilinear identity eq. ( A.1U| ). This 
gives: 



CO u Res W { k %(p', 

]T Res Wi%(p', p K )dS p >M (5.32) 



W^Pk) = - Res W^ x {p',p K )dS p ,, {p) 
p'^p 

r(h) 



We now introduce eq. ( |5.31| ) inside this formula and keep only terms which have 
poles at the branch points: 



wt%( P ,p K ) = 

- E h m 2 E s Res^ wl h - m \p')Rt\ P '; p K )dS^ {p) 

- E™=o EjZlm+rto S s Res p ,^ aa R { ™\p'\ Pj)W£^$(j/, p K -j)dS p > i0 (p) 

- J2 S Res p /_ as R£_^> (j/; p', p k )dSp,, (p) 



For convenience, let us note: 

Wi°\p)=0 (5.34) 
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Then, the recursive definition of W^ +1 (p,P^r) reads: 



Wffi(p,P*) = 

Ei=l Ern=0 Sj=o,m+j^O 

+ Efei E s Res p /^ as i^+i P'; Px)^p', (p) 

(5.35) 

5.2.2 A recursive formula for i2£ 

The second step consists in the derivation of an equivalent formula for Ft%. . We 
proceed in the same way as for the genus case: we use the rational properties of 
some of the correlation functions to write the recursive formula, with the aid of a 
relation similar to eq. (|5.22|) . 
Let Gt\x(q),y(r))he: 



Gt\x{q\y{r)) = {y{r) - y{q))Ui%, y(r); Pk) + ^ ^^'^ 

m=l j=0 U - L \'d) 
fc-1 j 

+ Z ~ch^q) Uj ( q ' PjWk-j+iiq, Pk-j) 

(5.36) 

The loop equation eq. ( |5.29| ) shows that G^\x(q), y(r)) is a rational function in 
x(q) and a polynomial in y(r). 
Thus, one has: 

Gf (stf), v(p')) = Gf\x(p),y(p i )) (5.37) 
which can be written: 



{vw)-y{p)m \p,y\p)',PK) = Em=oEL — : 



dx(p) 



dx 



We now establish a relation similar to eq. ( |5.22| ) in order to present our recursive 
formula in such a way that it can be graphically interpreted. 
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In order to achieve this aim, one has to determine an explicit intermediate for- 
mula for uji h \p,y;pK)- Let us assume that (for the proof, see Appendix C): 



Pk, 



E y (x,y(p 1 )) s -^min(d2,k+h) ^ sr^k+h v-< J_ 

y(p l )-y(p) ^r=\ l^Ki (J... [j K r =K 2^h a =0 l^k a = \K a \ i'jo^j'eWHi} Q (-539) 

wfel(P 1 .PJCl V 1 ' 1 ,• V 1 '* 1 " 1 ^ 1 ) (lla=2 Wfca+xO^ " >P*« J"' 1 ,• V"'^ ~ ^ 1 )) 

where Q = Y\ a {k a — \K a \)\ is a symmetry factor and one has the following 
constraints: 

• Y,o{K + k a ) = h + k; 

• < \K a \ < k a 

One should note that the only external parameter entering these constraints is 
k + h. 

It is now possible to derive an equality equivalent to eq. (|5.22|) . One shows - in 
Appendix D - that: 

E E w%}( P , Pj )u^ m) Mpy,PK-j) + u^( P ,y(p);p,Pk) 

m=0 j=0;mjj^kh 

= EE E w^l(p\ Pj )utr ) (p l Mpy,PK-j) 

i=l m=0 j=0;mj^kh 

+ Ilui h + - 1 1 \p\y(py,p\Pk) 

(5.40) 

This equality allows us to write: 



(j/(p*) - y(p))u { k h ' \p,y{p i )] 



^ h ^k ^ w^(p\ PJ )ui h - m) (p l >y(p^PK-j) . , 

+ X m + Wi%(p\ p K )E y ( X , y(f)) 



dx(p) 

That is to say: 



i%£- 1 >( P v, Pfc ) tOg%P£) 

Zj, ^ ! (y(p i )-y(p))dx(p) ~T~ (y(pi)-y(p))dx(p) 



(5.42) 
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5.3 Diagrammatic solution: a cubic theory 



This section is the principal part of the article. We define a correspondence between 
the correlation functions and a system of Feynman-like graphs. To every fc-point 
function of genus h, we associate a graph with k external legs and h loops and 
eq. ( |5.35|) and eq. (|5.42| ) become two relations describing these graphs as functions 
of graphs with less legs or loops thanks to some rules we introduce in this part. 

First of all, let us represent diagrammatically eq. ( |B.6|) and eq. (|B.16|) as the 
propagators of the theory: 



W 2 (p,q) = P 



(5.43) 



and 



(5.44) 



These two diagrams represent the basis of the whole representation: they allow 
to draw the k > 2 correlation functions. 

Note that the second propagator can also be seen has a vertex of valence 2, and 
this is the way it will be presented in the diagrammatic rules. 

Let us now introduce the whole diagrammatic representation: 



Let , and respectively, be represented as white and black disks with h 
holes and k external legs (remember that Wj^^ is the generating function of discrete 
surfaces with k + 1 boundaries and h holes): 



r(h) 




(5.45) 



Rl ih \p,PK) := 




(5.46) 



Let us introduce also the following propagators and vertices: 
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non-arrowed propagator: 



1 := W 2 (p,q) 



arrowed propagator: 



P ► q := dS q , (p) 



Residue cubic-vertex: 



colored cubic-vertices: 



nl 




(l-<5|, m )(l-<* m ,i)(l-<S M ) 
(y(p t )-y(p l ))dx(p) 



2-valent vertex: 



One can now simply interpret the recursion relations eq. (|5.35| ) and eq. ( 5.42|) in 
terms of diagrams. 

The relation eq. ( p.35| ) reads: 




d'2 h k 

EE E E p 

i=l m=0 j=0,m+j^0 JeKj 



d 2 



+E^ 



i=l 




(5.47) 



K-J 



And given lower order R]^ m ^s and w/ s, one can obtain diagrammatically 
by writing eq. ( |5.42| ): 
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(5.48) 



From these diagrammatic relations, one can see that W^ +1 is obtained by the 
summation over all diagrams with 1 root, k leaves and h loops following the rules: 

• The vertices have valence 2 or 3; there are 2h + k — 1 trivalent vertices; 

• The edges, are arrowed or not, the arrowed edges are waved or not; 

• The subgraph made of arrowed edges forms a skeleton tree (i.e. a tree whose 
vertices have valence up to 3); 

• from each trivalent vertex comes one waved and one non-waved propagator; 

• two vertices linked with a waved propagator have different indices; 

• the k leaves are non- arrowed propagators finishing atpj's (i.e. B(.,pj)); 

• the root is an arrowed non waved propagator starting from p. 

A practical way to draw these graphs is to draw every skeleton tree of arrows, 
put k non arrowed propagators as leaves, close it with h non arrowed propagators 
linking one vertex to one of its descendents in order to obtain h loops and then put 
waves so that from each trivalent vertex comes one waved and one non-waved arrow 
with the possibility that every waved arrow leads to a bivalent vertex. 

Remarks: 

• The order for computing the residues is following the arrows backwards from 
leaves to root. 

• Wk+i is symmetric in its k + 1 variables, although it is not obvious from this 
representation. 

• There is no symmetry factor arising in this representation unlike [13]. 
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5.4 Examples 

Let us briefly show some diagrams for small h and small k. 
5.4.1 Leading terms: tree level 

We begin by the leading terms of the first correlation functions, i.e. for h = 0. 
• k = 3: 




Ej=l Es ReSp'.^^ 



and 



+ p 



B(p'',pi)B(p',p 2 ) , B(p'\p 2 )B(p', Pl ) 
(y(p' l )-y( P '))dx(p') (y(p' i )-y(p'))<^(p') 



dS p/>0 (p) 



(5.49) 



^' (0 W,p 2 ) = E-Li 



<*2 




+ P' 



+ p 




(5.50) 



Let us show that W / 3°^(p,pi,p2) is indeed symmetric in pi, p 2 and p 3 . 

For every branch point a, let g be the only q % such that <ir(g) — > when g — > a. 



Wf W,P 2 ) = E^E S Res„ g(9 ' Pl) S^Sa g(9 ' pl) ^,o( P ) 



V d2 V Rp« Rp« B(r,p 1 )B{q,p2)+B{r,p 2 )B{q,p 1 ) ,g / \ 

Li=lLs ™^^ s ^r^?' (2/(r)-2/(?))(:r(r)-:E(g))<fa:(?) uo 9,°W 



= J2 S Res 
= Es Res 



B(r,p 1 )B(q,p2)+B(r,p2)B(q,p 1 ) , g / \ 
9 ^a s "^r^ (y(r)-y <?))(a;(r)-x(g))dx(< ? ) ao q,o{P) 



q^a s 



= -E s Res, 



Res, 

B(q,Pi)B(q,p 2 )'dS q ^(py 

(y(q)-y(q))dx(q) 

B(q, Pi)B(q,p 2 )dS q} q(p) 



(5.51) 



(y(9)-3/( g ))<ia;(g) 

— V Rpo B{q,p{)B{q,p2)B(q,p) 

— l^s ^^q^a s dx(q)dy{q) 



which is nothing but the formula found in [30] and is a way of writing Rauch's 
variational formula. 
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• k = 4: 



wi Wi,P2,p 3 ) = ££iE?LiP 



"T 2^i=l l~ij=\ P 



(5.52) 



"t" 2^i=i 2^j^i=i P 



+ ( permutations of {^1,^2,^3} ) 

One has to consider all the permutations on the external legs. Thus, is the 
sum over 18 different diagrams. 

5.4.2 Topological expansion: one and two loops level 

Consider now the first non planar examples beginning by the simplest one, the one 
loop correction to the one point function. 
• k — 1 and h — 1: 



W[ 1] (x(p))dx{p) 




d 2 




(5.53) 



=1 s 



yW) - y(q) 
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One can check that this is identical to the result of [17]. 
• k = 2 and h — 1: 




(5.54) 



Analytically, this reads: 

i=l Z^j6[l,d 2 ]-{i} Z^s aes p'^« s (j/(p' I )-S/(p'))(s/(p' i )-S/(P ,j ))^ 2 (p') 

[B{p', Pl )B{p'\p'i) + B(p' i ,p 1 )B(p',p'i) + B(p',p H )B(p 1 ,p^)} 

I v^(i 2 V Rpc p„„ dS p',o(P) 

tz,!=i Z.piLs,! J^-e&p'^a s -rxebp//_^ at ( 2/ ( p 'i)_j / (p'))(j / (p'o')_j / ( p »)) ( i x (p') ( i : i;(p») 

[s(p^p 1 )s(p^p^)^ >0 (p' ^ )+ J B(p^p 1 )s(p'^p^)^ I0 (p') 

• /c = 1 and h = 2: 
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(5.56) 



6. An effective non cubic theory 

The Feynman-like graphs described up to now correspond to cubic vertices only, but 
the price to pay is the introduction of auxiliary functions Ft^ h \ Nevertheless, in 
order to study some problems, this property is not needed and one may prefer an 
effective diagrammatic representation for only but vertices with valence up to 
d>2 — 1. This section is dedicated to building such a diagrammatic representation. It 
consists in resumming the linked waved vertices into one multivalent vertex: 




(6.1) 



6.1 Leading order: Genus 

We have already written the equations necessary to define this effective theory. Let 
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us consider eq. (|5.7|) and eq. (|5.16|) : 



w k+lip , PK) = - E Res g E w ™«> *MM 

(6.2) 



, E(x(p) : y)dx(p) ^ x ^ 

U k {p,y,PK) = — - — — E E 



dz 



W\ Kt \+l(p jt ,PK t ) 



y - yip) 

This second equation taken for y = y{p) reads: 



Y TT 



(6.3) 



Uk(p,y(p);i>K) 



W lKt \+l(p jt ,PK t ) 



d,2 d,2 r 

= y y y w 

E y (x(p),y{p))dx(p) ^ Kr^.JyK r =Kh+h^3T=\t=\ (vip) - y(P 3t )) dx iP) 



(6.4) 



Introduce it in eq. (|6.2|) and get a closed recursive formula for the W^s: 



Wk + i(p,PK) — ~ J2 S R es p'^a s ^t=lJ2K UK 1 U...UK r =K^l^j 2J L...jLj r =l 

W ]Ko \+i(p', Pk ) n t =i (^^$^f^dS P ',o{p) 



(6.5) 



Let us introduce the following Feynman rules: 



non-arrowed propagator: 


P 


q := W 2 {p,q) 


arrowed propagator: 


P 


► — q := dS q , (p) 


r+2 - vertex 




q 


(1 < r < d 2 ) 






with one marked 




\ ~~T q ]3 ^ t=1 {y(i)-y(<i H ))dx{q) 


edge: 










q J ' 



Remark that one leg of the multiple vertex is marked: on this leg, there is no 
summation over the different sheets. 



Using these rules, one can diagrammatically write the recursive relation as fol- 
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lows: 




d 2 

E E 

r=l K UKiU...UK r =K 




From this relation, one can see that Wk+i(p,PK) is obtained as the summation 
over all trees with k + 1 external legs and following the rules: 

• The vertices have valence r+2 such as 1 < r < min(k — 1, c?2); 



T/ie edges are arrowed; 



• One of the legs of each vertex is marked 

• The k leaves are non arrowed propagators ending at pj 's; 

• The root is an arrowed propagator starting from p. 



The drawbacks of these effective rules induced by the existence of multivalent 
vertices is balanced by the simplicity of the vertices and the absence of different 
propagators. 



6.2 Any genus h 

Let us now study the extension of this theory to any genus. 

Once again, the fundamental equations have already been written. Let us recall 
to mind eq. ( |5.35| ) and eq. ( |5.39| ): 

W { k %(p,p K ) = 

=0,rn+j^oY,s R-eS p /^ as Ey (x(p'),y(p')) 

) M-' j +Kv'^K-j)dS pl , a {jp) (6.7) 

_ p U ( k + l \p',y{p')\P'VK) iq I X 

l^s -n-eS p /^ as E y (x(p'),y(p')) aa p',a\P) 

and, for i ^ 0: 
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u ( u\pMv%vk) = 

E y (x,y(p 1 )) T—^min(d2,k+h) ^ sr^h sr^k+h J. 

y(p % )-y(p) 2 ^ r=1 ^Ki\J...\jK r =K2^h a =o2^k a =\K a \2^j a!0 ^j a f g ,ell,d 2 ]-{i}Q /qq\ 

vOi (pi ' p ^' pJ1,1 '---' p(3l ' fclH ^ 

In order to introduce this second formula inside the first one, one has to use the 
interpolation formula to consider the case where i — : 

t/ ; (m) (p,y(p);PL) 

E y (x(p),y(p)) 

i-^min(d2,l+rn) ^ ^ m ^l+m 1 

^r=l l^L 1 \J...\jL r =Ll^ ma =0l^l a ^L a \^ji^...^jr€[l,d2]n (6.9) 

rfx(p) r - i - 1+ ^' Q (s / (pJi,o)-j / (p))J] a(3 ( 3/ (pn,o ) „ y( y Q ,/3) 

Recursively, it is easy to check that it can be written: 

Ui m \p,y{p);pL) 

E y (x(p),y(p))dx(p) 

min(d 2 ,l+m) m l+ m ^ 

^ < rv 

r=l L 1 \J...\jL r =Lm a =Ql ol =\L a \j aiS ^j al0l e[lM 

iii da:(p)'--l^l+i nfc) LQl (l/(p) - 2/(p f -" ) ) 

(6.10) 

where Q' is some other symmetry factor depending only on the same parameters 
as Q. 

One is now able to write an explicit recursion formula for the W% s that can 
be graphically represented with the Feynman rules introduced in this section. The 
introduction of eq. ( |6.10| ) in eq. ( |6.7| ) gives: 



di h k+h i 

-ER es E E EE E 7^7 



wt+i(v' jafi iVK a iP' ja <\ • • • 



(6.11) 



E ^s ^(^(pO;^)^ (p) 

T^ a - Ey(x{p') } y{p')) 
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That is to say: 



P 




Remark that we have splitted the diagrams in the RHS in order to reproduce 
the recursion relation. Nevertheless, the first term in the RHS is nothing else but a 
particular case of the second term where the marked leg of the vertex is left alone 
inside one of the W's. 

Hence, the h-th order expansion term of the correlation function is obtained 
as the summation over all Feynman diagrams with k + 1 external legs and h loops 
following the same rules as exposed in the genus case, i.e.: 

• The vertices have valence r+2 such as 1 < r < d 2 ; 

• The edges are arrowed or not; 

• One of the legs of each vertex is marked; 

• The subgraph made of arrowed edges forms a skeleton tree; 

• The k leaves are non arrowed propagators ending at pj 's; 

• The root is an arrowed propagator starting from p; 

• a non arrowed edge links a vertex to one of its descendants along the tree. 
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6.3 Examples 

Let us review some simple examples of this description. 
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7. The gaussian case: the 1-matrix model limit. 

In this section, we are interested in the special case where d 2 — 1, i.e. one has a 
gaussian potential in M 2 . This situation is very important because it links our results 
to the 1-matrix model studied in [13]. Indeed, when one of the potentials is gaussian 
- V 2 for example -, the integration over one of the variables - M 2 in this case - is 
gaussian and can be straightforwardly performed without giving any contribution 
to the formal expansion. Then, the 2-matrix model with one gaussian potential 
V 2 (y) — ^U 2 is equivalent to the 1-matrix model with a potential V — V± — J^-. We 
check in this part that our results coincide with the ones obtained directly from the 
1-matrix model in [13]. Actually, it is a good way to better understand the structure 
obtained. 

In this case, the Riemann surface is an hyperelliptical surface with only two x- 
sheets. The equation x{p) = x has only two solutions. Let us call them p and p, i.e. 
p° = p and p 1 = p. They obey the following relations: 

x(p) = x(p)y(p) = -y(p) (7.1) 

The algebraic equation generating the Riemann surface reads: 

E{x(p),y(r)) = -g 2 {y(r) - y(p))(y(r) - y(p)) = -g 2 (y(r) 2 - yipf) (7.2) 



One can also remark that: 



Ukip, V, Pk) = g 2 W k+1 {p, p K ) (7.3) 

That is to say: 

no ( \ Uk(p,y(p);i>K) w fc+ i(p, p K ) (7A] 

k[P ' PK> E y (x(p),y(p))dx(p) ~ " 2y(p)dx(p) { ' } 



So that: 



R k (p,PK) = R k ( P , PK )= 2y{p)dx{p) (7.5) 



Diagrammatic rules. 

One can now study how the diagrammatic rules introduced earlier behave in this 
limit. 

• The cubic rules 

Because V 2 is gaussian, the Feynman rules become: 



33 



non-arrowed propagator: 



arrowed propagator: 



Residue cubic-vertex: 



simple vertex: 



q := W 2 (p,q) 



q := dS q!0 (p) 



■ X/s R&Sq—Hi B 



-p : = 



2y(p)dx(p) 



The last component of the Feynman diagrams, the colored cubic-vertex, im- 
plies three different x-sheets. Because there exists only two such sheets in the 
gaussian case, this vertex vanishes: 

. m 



(7.6) 



Considered that the bivalent and trivalent vertices only appear together, one 
can merge them into one whose value is equal to —J2 S R- es g->a s 2 y {q)dx{ q ) ' anc ^ 
one recovers [13]: 





9 




(7.7) 



The effective theory 

The effect of the gaussian limit on the effective theory is to make it cubic. One 
obtains the following rules: 



non-arrowed propagator: 



q := W 2 (p,q) 



arrowed propagator: 



q := ds q M 



cubic vertex 
(only for r=l): 




» 2y(q)dx(q) 
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Hence, the two theories turn into only one cubic theory in this limit which is the 
one derived in [13]. Indeed, the corresponding recursive relation appears to be: 



<W) = - E Res W ^f^M 

h j= okm^i^ 1 2y(q)dx(q) 

Remark: 

Diagrammatically, this limit can be easily interpreted. Starting from the general 
cubic theory, in order, to obtain the 1-matrix model graphs from the 2-matrix model 
ones, one only has to take the length of the waved propagators to 0. In this case, 
the graphs containing at least one colored vertex vanish. 

Everything works as if the waved propagators of the 2-matrix model were un- 
stable particles which decay into stable ones represented by non-waved propagators. 
Then the 1-matrix limit is obtained by taking the life time of these particles to 0. 



One shall also note that there is no symmetry factor in the 2-matrix model graphs 
of the cubic theory whereas there are not well understood ones in the 1-matrix case. 
The derivation of the 1-matrix model as a limit exhibits how these factors arise. 
They come from the same contribution given by different diagrams in this limit. 
This observation exhibits how the 2-matrix model seems more fundamental. 
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8. Conclusion 



In this article, we have generalized the diagrammatic technique of [13] to compute all 
non-mixed correlation functions of the 2-matrix model, to all orders in the topological 
expansion. 

The result can be represented diagrammatically, with some cubic Feynman rules, 
which are just convenient notations for writing residues on an algebraic curve and it 
is not clear whether there exists a field theory giving rise to these graphs or not. 

This shows that the method discovered in [13] is very universal, i.e. it works for 
all algebraic curves, not only hyper elliptical curves. 

The future prospects of that work are to find the diagrammatic rules for com- 
puting the free energy to all order in the topological expansion, and also all mixed 
correlation functions (using the result of [14]). Another possible extension is to work 
out the multimatrix model, i.e. the chain of matrices as in [20], and in particular the 
limit of matrix quantum mechanics. We believe that this technique could apply to 
many other integrable models. 

Another question, is to understand the limit of critical points, i.e. when some 
branch points and double points start to coalesce. It seems that the diagrammatic 
technique should just reduce to consider only residues at branch points which become 
critical. One may expect to recover some relation with the Kontsevich integral, in 
relationship with KP integrable hierarchies. 
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Appendix A Needed tools of algebraic geometry 



We review here some definitions and properties all along this article. 

Behaviors at oo. We see from eq. (|2.17|) , that at large x, we have y ~ V[{x) — 
^+0(l/x 2 ) in the x-physical sheet, (resp. at large y, we have x ~ V^y) — ^+0(l/y 2 ) 
in the ^/-physical sheet). This means that the functions x(p) and y(p) have two poles, 
oo + and oo_ on £. The function x(p) has a simple pole at oo + and a pole of degree 
d 2 at oo_, while the function y{p) has a simple pole at oo_ and a pole of degree di 
at oo + . We have: 

V® P -Z> + V ^ x ^ -^pj + °(V^(P) 2 ) (A.l) 
*(P) ™ V MP)) + O^/Vipf) ( A - 2 ) 



p- 



In particular: 



Res ydx= Res xdy = 1 (A. 3) 



Genus and cycles. The curve £ is a compact Riemann surface with a finite 
genus g < d\d% — 1. If g — 0, £ is simply connected, and if <? ^ 0, there exist 2g 
linearly independent irreducible cycles on £ , such that by removing those 2g cycles 
we get a simply connected domain. It is possible to choose canonically the 2g cycles 
as Ai-, &h i — 1) ■ • • j SS such that: 

Ani j = fl , Bi n Bj = o , An ^ = <j y (A.4) 



Branch points. The x-branch points dj, z = 1, . . . , di + 1 + 2g, are the zeroes 
of the differential dx, respectively, the ^/-branch points bi, % — 1, . . . , d\ + 1 + 2g, are 
the zeroes of dy. We assume here, that all branch points are simple and distinct, i.e. 
that the potentials are not critical. Notice also, that E y (x(p),y(p)) vanishes (simple 
zeroes) at the branch points (it vanishes in other points too). 

Bergmann kernel. On the Riemann surface £, there exists a unique Abelian 
bilinear differential B(p, q), with one double pole at p = q, such that: 

B(p,q) ~ d <P) dx ^) +finite and v , I B (p,q) = (A.5) 
p^i (x(p) -x(q)) 2 JpeA, 

It is symmetric: 

B(p,q) = B(q,p) (A.6) 

Its expression in terms of theta- functions can be found in [24,25], it depends only on 
the complex structure of £. 
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Abelian differential of third kind. 

On the Riemann surface £, there exists a unique abelian differential of the third 
kind dS q:r (p), with two simple poles at p = q and at p — r, such that: 

Res dS q>r (p) = 1 = - Res dS q:T (p) and Vi jf dS q:T (p) = (A.7) 

We have: 

dS,, r (p)= f S(p,g') (A.8) 

J g'=r 

where the integration path does not intersect any Ai or £>j. 

dS q>r (p) is a differential on £ in terms of p, but it is a multivalued function of 
q (and of r). After crossing a cycle Bi, it has no discontinuity, and after crossing a 
cycle Ai, it has a discontinuity: 

disc (dS qtr (p)) = dS q r (p) - dS q _ >r {p) = f B(p, q') (A.9) 

Jq'eBi 

Note that the discontinuity is independent of q. 
Riemann bilinear identity. 

If u; is a differential form on £, such that § qe ^ uj{q) = 0, we have: 

9 r 

H Res u(q)dS q , r (p) =Y,f disc A(^(9)^ 9 ,r(p)) 

i q Zl i=l Jq£Ai 

-H f „ discB^u^dSg^p)) 

9 y 

= 12? ^iq) disc A ^dSg tr (p)) 
i=l JqeAi 

9 j- 

= ^disc A (c£V(p)) f ^(Q) 
i=i JqeAi 

= (A.fO) 

where the LHS is the sum over all residues on a fundamental domain, the poles Zi are 
all the poles of uo as well as the pole at q — p. This identity is obtained by moving the 
integration contours on the surface, and taking carefully into account discontinuities 
along the nontrivial cycles (see [24,25]). 

Appendix B Two points function in the planar limit 

We present here a new derivation leading term of the 2-point function's leading term 
W 2 (pi,p 2 )- 

This case is of special interest because it represents some initial condition for the 
diagrammatic rules. In fact, the two correlation functions W 2 (pi,P2) and Ui(pi, y; p 2 ), 
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are the basis of the whole structure of the wj^'s. Moreover, it allows us to show 
through a simple example the way we proceed further for the general case. 

We first rederive the well known result that the two point function is nothing 
else but the Bergmann Kernel (see [5] for instance). 



Let o G £ be an arbitrary point on the Riemann surface. Since the Abelian 

dx(p) 
x(p)-x(q) 



differential of the 3rd kind defined in eq. ( |3f ) dS q:0 (p) behaves as x ^^l q \ when 



q — > p, one can write the Cauchy formula under the form: 

W 2 (p, Pl ) = - Res dS q , (p)W 2 (q;pi) (B.l) 



q->p 



One can see from eq. (|5.2|) with k — 1, and from eq. (|4.17|) , that the integrand 



in the RHS has poles only for q — > p and q — > pi, Since W 2 has vanishing ^4-cycles 
due to eq. ( |4.17| ), we can use the Riemann bilinear identity eq. ( |A.10| ), and get: 



W 2 {p, P i) = Res dS q>0 (p)W 2 (q;pi) (B.2) 



q^pi 



For k — 1, eq. ( |5.2| ) reads: 

E y( x (p),y(p)) W 2 (p,Pi) = -Pi(x(p),y(p);p 1 )dx(p) 

, ( U (pi,y(p)) dx(p) \ 

~rCLp 1 



(B.3) 

and thus we have: 



x(p) — x(pi) dx(pi) J 



W 2 (p;px) = R£sdS qt0 (p)W 2 (q;pi) 

™ E y (x(q),y(q)) 

+ Res dS 9 M P ^*(f)-*M ^)J 

™ E y (x(q),y(q)) 

(BA) 

Since Pi(x(q),y(q);pi) is a polynomial in and it has no pole at q — p\. 
For the second term we use eq. 



W 2 {p\pi) = Res dS q>0 (p)W 2 (q;pi) 

= d Res dS (p) E(x(p l ),y(q))dx(q) 

P1 n g '° (x(q)-x(p l ))(y(q)-y(p 1 ))E y (x(q),y(q)) 
= d pi dSp lfi (p) 
= B(px,p) 
(B.5) 
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We thus recover the well known result: the two-points function is equal to the 
Bergmann kernel on the Riemann surface corresponding to the algebraic equation 
E(x,y) = (cf [5,20,27,30]). 



(B.6) 



Let us now compute Ui(p,y;pi). For k — 1, eq. (|5.1| ) reads: 

(yjr) -y(q))Ui(q,y(r);pi) _ W 2 (q;pi)U (q,y(r)) pff s / \ x 
dx(q) ~ dx{qf nWhVVhPi) 

Pl \(x(q) - x(pi)) dx(px) J 

take it for q = r = p % : 

W 2 {p i ;pi)U {p i ,y{p i )) i i 

= ^-p Pi(x(p),i/(p);pi) 

+d f Ua&uyjp*)) \ (Bg) 



using that = x(p l ), we have: 



= ^ Pi{x(p),v(pY,Pi) 

+a ^\{x{p)-x{p 1 ))dx{p l )) ^ 



Now, write eq. (|B.7|) with q = p and r = p % : 

(y(jt)-y(p))Ui(p,y(jt)]Pi) W 2 (p;pi)U (p,y(p i )) 



dx(p) dx{p) 



2 



- Pi(x(p),yip' l );pi) 



+d i f Uojpuyjf)) \ , B 1Q) 

Pl \(a;(p) -x(pi))da;(pi)y 

and insert eq. ( p.9|) , you get: 

W 2 (p l ; Pl )U (p\y(p 1 )) 



(y(p l ) - y(p))Ui(p, y(p l );pi) 



dx(p) 
W 2 (p; Pi)U Q {p,y{p i )) 



dx(p) 

(2J.11) 

Using eq. flOft , i.e. U (p,y) = ^^ dx(p), this implies: 

(y(p l )-y(p))Ui(p,y(p i y, Pl ) = W 2 (p l ;p 1 )E y (x(p l \y(p i )) (B.12) 
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Since Ui(p, y;pi) is a polynomial of degree <i 2 — 1 in y, we can reconstruct it through 
the interpolation formula: 

E(x(p),y) A 1 (yip 1 ) -y(p))U 1 (p,y(p i );pi) 



(5.13) 



i.e. 



Ejxjp),y) ^ W 2 jp\ Pl ) 

(y - yip)) h» _ 2/(p ) 

and in particular, at y — yip), we have: 

R^p^dxip) = = ^ — — B.15 

Ey{x{p),y(p)) fzi yip) - yip 1 ) 

and for i ^ 0, we have: 

R^pxjdxip) = = — — — B.16 

E y {x(p), yip 1 )) iyip 1 ) - yip)) 

Appendix C Computation of Eq. (5.39) 



In this appendix one proves recursively eq. Q5.3S ) for any k and h. 

Let us suppose that this formula is know for any C7j with m < h — 1 and for 
any Uf 1 with I < k — 1. One writes it: 



C// m) (p,2/(p l ');Pi) = 

,k+h) sr-\ \r^m v-v£+m v-^ J_ 



y{p l )-y{p) ^ r=1 ^L 1 \J...[jL r =L 2^m a =o 2^i a= \L a \ 2^j a> p^j al /3 ieli,d2}~{i} n 

Let us introduce some shortened notations so that one can write this proof in a 
few pages. 

Considering the sum on the RHS of eq. ( |C.1| ), one can see that there are two 
different kinds of terms: 

• If lx — \Li\, one can factorise the term W^^ip 1 , PlJ. Let us note the sum of 

these terms W(p l ,Pl)W(pl,p') where we have noted W instead of W^ 1 ^ to 
indicate that these are formal notations; 

• the other terms correspond to the sum over all l\ ^ \L±\. Let us denote them 
by W(p i ,pL,P j )W(p L ,P j )- 
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Using these notations, one can shortly write eq. ( |C.1|) : 

ut\p,y(p i y,PL) = w(p\ PL )w( PL ,p i ) + w( P \p L ,p j )w(p L ,p j ) (C.2) 

Thus the interpolation formula gives: 

ul m \p\y(p%p L ) = W(p L)T P) + W(p,p L )W(p L ,pi) + W(p,p L ,pi)W(p L ,pi) (C.3) 

where the first term corresponds to the sum where all jps are different from i 
and and there is no whose argument is p or p l . 
On the other hand, one knows the relation |5.38| : 



m=0 £^j=0 (yip^-yip^dxip) 

_s r h w j7i (p-p j) u i h -r ] (p>y(p i y<PK-j) (C4) 

2^ m =0 2^j=0 (j/(p 4 )-J/(p))rfx'(p) V ' 1 

(l/(p*)-S/(p))da: (l/(p')-«(p))<te 

Remark that the terms in the RHS of this equation correspond to the criterion 
of the hypothesis and one can then express them as a product of W's following the 
notations introduced earlier. This reads: 

ui (p,v(j?)]pk) = 

W(p\p K )W{p K ,p>) + W{p i ,p K )W{p,p K )W{p K ,p j ) 
+W(p i ,pK)W(p,pK,P j )W(p K ,p j ) - W(p,p K )W(p\pK)W(p K ,p j ) 
-W{p,pK)W{p\pK,P j )W{p K ,P j ) + W(p i ,p K ,p j )W(p K ,p i ) 
+W(p\p,p K )W(p K ,p j ) + W(p,p K )W(p\p K ,p j )W(p K ,p i ) (C.5) 
+W(p\p,p K ,p j )W(p K ,p j ) + W(p,p K ,p>)W(p i ,p K ,p j )W(p K ,p j ) 
-W(p,p\p K )W(p K ,p j ) - W(p i ,p K )W(p,p K ,p j )W(p K ,p j ) 
-W{p,p i ,p K ,p j )W{p K ,p i ) - W{p i ,p K ,p j )W{p,p K ,p j )W{p K ,p j ) 
= W{p\p K )W{p K ,p>) + W{p\p K) p j )W{p K) p>) 

So one has proven the formula for llj^K 

Because this formula is true for h=0, it is true for any k and h. 



42 



Appendix D Derivation of Eq. (5.40) 

One wants to show that: 



Eti eLo zU;mm h w$\ { P \ p.j)utT ] tf, v(p);pk-j) 



u[ h , 1) (p,y(p);p,Pk) 



dx 



(D.l) 



+ Y? 2 ^ S±L 



(p\y(p);p\Pk) 



dx 



Let us compute the difference D between the two sides of the equation by the 
introduction of eq. ( |5.39| ) written with some few different notations which are defined 
as follows: 

• I = r + h -J2 a h a ] 

• u p = TLi(k e -\K e \)-/3. 



One can then write: 
D 



m=0 ^j=0;mjy^kh 



W%}(p,pj)E y (x,y(p))j:« 



d> 



x Z^r=l l^o=l 2^h o =0 2^k 



1 y(p)-y{p l ) 

h-m sr^k+h-j-m ^ 1 
=0 2^j 2 ^...^j l e[l,d2}-{i} 2^K 1 \J...[jK r =K Q 



I'V, 



X- 



Oi) 
fci+i 



Ed 2 ^-^k 
i=l l~,m=0 2~ij=0 



Y\\ =2 y{p 1 )-y{p 3 -<) 



7=2 ■ 
r (m) 



v v^^2 v^r ^h—m •^k+h—j—m ^ y^ i 

x 2^r=i 2^o=i 2^h D =o 2^fc D =o 2^j 2 ^...^j ; g[i,d 2 ]-{j} 2^ (J... (J # r =K- ^ 



X 



<|\ (PiPK 1 ,p>r + l ,...,p>r +U1 ) ( n ^ a frfr ^ .p'^-l 



nU2^(p)-y(p j7 ) 



+-g y (g,y(p))Eiii 

x 2^r=i 2^o=i 2^h o =0 2^k =o 2^j 2 ^...^j t e[i,d 2 ]-{i} luKi \J...\jK r =K n 



x 



n 7 =2^(p i )-^^ 7) 



W^Vi^'P^i ,P ir+1 ,-y r+tt i )W^ 1 (p,p>2 , PK2 ^+«i+i ) .„yr+« a ) 



_ E vi x ,y{p)) 



U' 2 y(p')-y^) 



A=\ y(p)-y(pi) 

ir^d, 2 ir^r v-*' 1- 1 v-~ > ' c +' 1 Y^ Y^ 1 

x 2^r=i 2^o=i 2^h o =0 2^k o =0 2^j 2 ^...^j l e[i,d 2 ]-{i} z^Ki \J... \J k t =k n 



x 



w < k h ^ 1 {p\p,PK 1 ,P 3r +\ 



/9=2 " fcg+1 



(p^,Px„p Jr+t '' 3 - 1+1 ,..,p Jr+t '' 3 ) 



n'= 2 2/(p)-?/(p j7 ) 



x n^ =3 w^^.p^y^"-^ 1 , • • • y^) 

The difference between the two first terms leaves only the terms corresponding 
to Mi 7^ in the first one minus U\ 7^ in the second one. 



U\= 2 y(p i )-y(p j ' 1 ) 



(D.2) 
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The difference between two last terms will allow us to compensate the preceding 
ones. Indeed, the terms with p l and p together in the same correlation function 
straightforwardly vanish and one gets the exact opposite to the two first terms re- 
maining. 

Thus D=0 and the equality |D.1| is proven. 
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